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ABSTRACT 

We present the spatially resolved emission line ratio properties of a ^10 10 M Q star-forming galaxy 
at redshift z = 1.03. This galaxy is gravitationally lensed as a triple-image giant arc behind the 
massive lensing cluster Abell 2667. The main image of the galaxy has magnification factors of 14±2.1 
in flux and ~ 2 x 7 in area, yielding an intrinsic spatial resolution of 115—405 pc after AO correction 
with OSIRIS at KECK II. The HST morphology shows a clumpy structure and the Ha kinematics 
indicates a large velocity dispersion with V max sm.(i)/a ~ 0.73, consistent with high redshift disk 
galaxies of similar masses. From the [Nil] /Ha line ratios, we find that the central 350 parsec of the 
galaxy is dominated by star formation. The [Nil] /Ha line ratios are higher in the outer-disk than 
in the central regions. Most noticeably, we find a blue-shifted region of strong [Nil] /Ha emission in 
the outer disk. Applying our recent H II region and slow-shock models, we propose that this elevated 
[Nil] /Ha ratio region is contaminated by a significant fraction of shock excitation due to galactic 
outflows. Our analysis suggests that shocked regions may mimic flat or inverted metallicity gradients 
at high redshift. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift — gravita- 
tional lensing: strong — ISM: jet and outflows 



1. INTRODUCTION 

One of the most challenging aspects of current galaxy 
formation and evolution theory is the modeling of feed- 
back processes associated with star- formation and/or 
nuclear activities. Feedback is considered to be the 
key solution to many observational puzzles such as 
the quenching of star - formation in ma ssive galaxies 
(e.g., (Bower et all 120061 : iDekel et alJl2009D. the origin of 



the mass-metallicity (MZ) re lation (e.g., iTremonti et al.l 
12004 iFinlator fc Pavel 120081) . the enrichment of the in- 
terstellar and interg alactic medium (ISM, IGM) (e.g., 
IVeilleux et al.l I2005L and references therein), the flat- 
ter central dark matter pro file than that predicted 
in ACDM simulation s (e.g., iFlores fc Primackl 119941 : 
iGovernato et al.l [20121) . and the problem of the forma- 
tion of bulge-less gala xies (e.g.. IGovernato et alj|2010t 
iPilkington et al.ll2011l ). 

Galactic-scale winds and outflows are observed in 
a wide range of galaxy types that are undergo- 
ing significant star formation. Galactic-scale out- 
flows have been frequ ently observed in (1) loca l star- 
burst dwarf galaxies (IMartinl 119991: iKirbv et al.l 120111 : 
Ivan der Wei et all 1201 ll) . (2 ) local ultra- l umin o us in- 
frared galaxi es (ULIRGs) (jRupke et al.l l200l 120051 : 
IMartinl l2005f ) , and (3) high-redshif t sub- millimeter 
galax ies (SMGs) (|Nesvadba et alj|2007t lAlexander et al.l 
I2010T) . Galact ic winds are prevalent at intermediate and 
high redsh i ft dS haplcv et al. 2003; Nesva dba et al.l 120081: 
ISato et al.1 120091 : IWeiner et al.l 120091 : iMartin fc Bouchel 
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12003 IMartin et al.ll20H iRubin et al.ll2010bL l20ll . 

In cosmological hydrodynamical simulations, imple- 
mentation of a constant wind speed and outflow mass 
fails to reproduc e the galaxy stellar mass functions 
(e.g., ICrain et al.l l2009f ) . A momentum-driven wind 
feedback model is currently favored by observations 
and can reproduce the stellar mass fun ction and the 
observed large-scale metal enrichment dMurrav et al.l 
20051: iQppenheimer k Da ve 2008; Finl ator fc Davel2008t 
Dave et al.l 1201 if h In the momentum-drive wind 
model, the efficiency of the winds or mass-loading 
factor (M w ind/ SFR) scales with the declining power 
of the mass or circular velocity (VjT 1 ) of the galaxy 
(M w i n d/SFR oc V~ r ). However, whether this simple 
scaling relation applies to individual galaxies of a wide 
mass range at different redshifts remains to be tested. 
The physical picture of how exactly the winds interact 
with and recycle the energy and chemic al elements to 
the IS M and IGM is still lacking. Indeed, Hopkins et al. 
(2012) have recently developed a more realistic stel- 
lar feedback model that includes a complex multi-phase 
structure that depends on the interaction between mul- 
tiple feedback mechanisms operating on different spatial 
and temporal scales. 

Note that past studies of gas inflows/outflows have 
mostly been based on slit-spectra of a small number o f 
individual galaxies (|Sato et alJl2009URubin et alj|2010al) 
or average com posite spectra from stacked low signa l- 
to- noise spectra (jWeiner et al. 2009; Rubin et al]l2010b| ). 
Recent studies using KECK /LRIS data have significantly 
expanded the sample of individual galaxies and prevalent 
outflows have bee n confirmed for hundreds of individua l 
galaxies at z ~ 1 (jRubin et al.ll2012t IMartin et al.ll2012l ). 

To provide a more concrete observatio nal baseline for 
the so phisticated feedback models such as Hopkins et al. 
(2012), the next leap forward is to look for wind sig- 
natures in a spatially resolved manner. Integral field 
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spectroscopic (IFS) observations are needed to sepa- 
rate t he various spatial components of the multi-p hased 
ISM (lArribas et al.l l200ll iRich et al.l 120101 [20T1 . Re- 
cent surveys of wide field of view IFS on local galaxies 
have dramatically expanded the quantitative analysis for 
galactic- wide winds and shocks dRich et al.l 120101 [20111 : 
lArribas et aI1l2008t ISharp fc BlancTH awthornl l2010l ). 

IFS studies have found widespread shock-excitation as- 
soci ated with gal actic outflows in massive local galax- 
ies (iFarage et all 120101 : ISharp fc Bland-HawthoTnl 120101: 
IRich et al I l201Ct 120111 ). Shocks are common products 
of winds and can also be induced by gas accretion and 
tidal flows in mergers (IVeilleux et al.ll2005t iFarage et al.l 
120101 IGenzel et al.l 120111: iSoto et al.ll2oH . Shock exci- 
tation contributes to the ionized state of the hot gas in 
galaxies, which also depends intimately on the chemi- 
cal abundance of the ISM. To identify the location and 
fraction of shock excitation across an entire galaxy, it is 
crucial to disentangle the shock induced ionization and 
energy dissipation from other mechanisms. 

Narrow line ratio maps can be u sed to probe 
the i onization nature of the sources (jVeilleux et al.l 
2002; IVeilleux & Rupkel l200l . Strong emission line 
ratio s such as the BPT diagr a m ([Baldwin et alj 
198lt IVeilleux fc qsterbrockl [l987t lAllen et all 119981: 
Kewlev et al.ll2001bl I2006T ) remains one of the most prac- 
tical methods to probe the ionization mechanisms at 
both low and high redshift. Combined with the most 
recent photoionization and shock models, strong line ra- 
tios such as [NiiJ/Ha, [Oiii]/H/3, and [Sn]/Ha offer a 
powerful tool to identify the excitati on mechan is ms in 
non-obscured regions. For example, IRich et all (|2010[ ) 
showed that in NGC 839, the [Nil] /Ha, [Om]/H/3, and 
[S 11] /Ha ratios are systematically lower near the nucleus, 
yet high along an outflow bi-cone. The region with high 
line ratios is remarkably well fitted by a shock model with 
slow velocities typical of galactic winds. 

At high redshift, the star-formation rate (SFR) of 
galaxies is su bstantially larger th an local star- forming 
(SF) galaxies (iNoeske et"aHl2007allbl : IWuvts et al.ll20TlE 
ISobral et al.l l2012bP). and galactic outflows are preva- 
lent and powerful flErb et al.ll2006l; ISwinbank et al.ll2007l: 
Weiner et all 120 09: Ru bin et al.ll2010bt lAlexander et al l 
2010t IGenzel et al.l 1201 If) . It is therefore of paramount 



importance to observe the masses and energies of high 
redshift winds to constrain the relations between the 
wind efficiency, SFR and galaxy masses on a broader 
scale. 

Because of the observational challenges, outflows and 
super winds at high redshift have mostly been reported 
for LBGs, SMGs, or in quasar absorption line sys- 
tems. Spatially resolved spectroscopy for high-z galaxies 
hav e been limited to the very massive st ar burst galax- 
ies (jCresci et al.ll2010l: IGenzel et al.lj20Tlt ). With seeing- 
limited observations, iCresci et al.l (|2010l ) reported three 
z ~ 3 UV-selected galaxies with inverted metallicity gra- 
dients (i.e., lower metallicity in the nucleus than in the 
outer-disk). The inverted gradient was explained by pri- 
mordial cold gas being funneled to the nucleus. These 
studies demonstrate the power of IFU studies in reveal- 
ing the underlying physical processes of the gas accretion 
at high redshifts. 

With the aid of magnification from gravitational lens- 
ing and adaptive optics (AO), we can begin to probe the 



detail s of fainter and less massive galaxies at high red- 
shift ([Nesvadba et alir2008HSwinbanldl2007tlJones et al.l 
12010b!: iSharon et all 12011 IWuvts et al.l 12012ft . The 
highly magnified multiple images in massive clusters can 
be magnified by more than ten, offering a rare opportu- 
nity to obtain spatially resolved observations. This tech- 
nique has enabled the first high resolution m easurements 
of metallicity gradients at high redshift (jStark et al.l 
2008; Uones et al.ll2010at lYuan et al.ll20TTI) . 

To determine the metallicity gradient evolutions with 
redshift, we are undertaking a metallicity gradient sur- 
vey on lensed galaxies behind galaxy clusters. Here we 
report observations of the z = 1.03 arc behind the mas- 
sive cluster Abell 2667. We find anomalous spatial line 
ratio variations indicative of large scale shocks. 

Throughout this paper we use a ACDM cosmology 
with H a = 70 km s" 1 Mpc" 1 , fi M =0.30 and fi A =0.70. 
At z=1.0, 1 arcsec corresponds to 8.1 kpc and a look- 
back time of 9.3 Gigayear (Gyrs). We use solar oxyge n 
abundance 12 + log(O/H) =8.66 (|Asplund et al.ll2005[ ). 

2. OBSERVATIONS, DATA REDUCTION AND 
ANALYSIS 

2.1. Observations 
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Fig. 1. — Three-color image of the A2667 giant arc at z = 1.03 
with FSTF814W (red), F606W (green), F450W (blue) filters. The 
white box on the bottom left shows the FOV and PA (37°) of the 
OSIRIS Jbb filter. The triple images composed of arc A, B, C. 
The IFU was centered on arc A (02000, feooo) = (23:51:40.023, - 
26:04:51.66). North is up, and East is to the left. 



The A2667 giant arc is a three-image lensing system. 
ICovone et al.l (2006) provided the spectroscopic redshift 
of z = 1.033 for the arc and the lensing models of 
the Abell 2667 cluster at z = 0.233. We observed the 
Abell 2667 giant arc with the OH -Suppressing Infra- R ed 
Imaging Spectrograph (OSIRIS; lLarkin et al.l I2006D in 
natural guide star Adaptive Optics (NGSAO) mode on 
the KECK II telescope at (a 20 oo, 5 2 ooo) = (23:51:40.023, 
-26:04:51.66). The observation was conducted on the 



night of September 20-21, 2011, using natural guide star 
(NGS) AO correction. The R=14 mag natural guide 
star was on the edge of the NGS ficld-of-view (FOV) 
of OSIRIS, restricting the positional angle (PA) to 37°. 
The 100 mas pixel scale was used. Since only broad band 
filters were available on OSIRIS at the time, we used the 
Jbb filter (spectral coverage: 1.18 to 1.416/im, spectral 
resolution: R ~ 3600, FOV 1.6" x 6.4"), which covers 
the Ha, [N n], and [S n] lines. The observations were con- 
ducted in a standard ABBA dithering sequence, with the 
A frame centered on the target and the B frame centered 
on an object-free sky region ^20 arcseconds North-East 
of the target. The same position angle was used for A 
and B (sky) frames. A total of 5 exposures were obtained 
for the target, with 900 seconds on each individual frame. 
The net on-target exposure time is 75 minutes. Figure Q] 
shows the OSIRIS FOV and locations of the A2667 arc 
on the HST F814W, F606W, F450W color-coded image. 

Before the science target exposure, short exposures of 
the tip-tilt star (TT) were taken to center the IFU. An 
offset was applied from the TT star to the science target. 
The offset was determined from astrometry corrected 
HST images with - 0.25" (i.e., ~ 2-3 OSIRIS pixels) 
accuracy. Gaussian fitting to the point spread function 
of the TT star yields an average FWHM of 0.13". The 
average FWHM of the sky OH lines is measured to be ~ 
2.5 A for our data and we use this value as the instru- 
mental resolution. In the analysis below, we deconvolve 
the line widths with this instrument resolution by sub- 
tracting the instrumental resolution in quadrature from 
the best-fitting Gaussian a. 

2.2. Data Reduction 

Individual exposures were fi rst reduced using the 
OSIRIS data reduction pipeline (jLarkin et al. 2006 ]). W e 
used the sky subtraction IDL code of iDaviea (|2007j ). 
which we have modified to optimize the sky subtrac- 
tion specifically around the wavelength region of Ha and 
[Nil]. The [Sir] lines are undetected due to strong tel- 
luric absorption. Since we are mainly concerned with the 
emission lines, a first order polynomial function was fit 
to the continuum for each spatial sample pixel (spaxel) 
and then the continuum was subtracted from each spaxel. 
The subtraction of the continuum helps to improve the 
removal of lenslet to lenslet variations. The final dat- 
acube was constructed by aligning the sub-exposures 
with the centroid of the Ha images and combining using 
a 3ct mean clip to reject cosmic rays and bad lenslets. Tel- 
luric correction and flux calibration were performed by 
averaging three AO standard stars (HIP220, HIP1272, 
HIP14719) observed immediately after the science ex- 
posure. Due to the low declination of this target, the 
airmass varied between 1.4 to 1.8 during the observa- 
tion. The uncertainty in absolute flux calibration is esti- 
mated to be ~20— 30%. However, the relative flux ratio 
of [Nil] /Ha is accurate to 10%. 

2.3. Analysis 

Our ana l ysis fo llows the procedures outlined in 
IJones et"aT] (|2010H) . Gaussian profiles were fitted simul- 
taneously to the three emission lines: [Nil] A6548, 6583 
and Ha. The centroid and velocity width of [Nil] A6548, 
6583 lines were constrained by the velocity width of 
HaA6563, and the ratio of [Nn]A6548 and [Nn]A6583 
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Fig. 2. — Observed spectrum for A2667 arc A, in the vicinity of 
Ha and [Nil] lines. This is an integrated spectrum from summing 
all IFU pixels with Ho detections above 5cr. The flux density is in 
units of 10 — 16 ergs s — 1 cm -2 A -1 . The raw spectrum is shown in 
black. Gaussian fits to the 3 lines are shown in blue. Note that the 
weaker [Nil] A6548 line is fixed to be the theoretical value (1/3) of 
the [N II] A6584 line. The velocity width of [N II] is constrained to 
be the same as Ha. Residuals are plotted in green. The locations 
of the strong sky OH emission (scaled and extracted from a pure 
sky frame) are shown in red at the bottom of the spectrum. The 
atmosphere transmission (black dotted line) is shown at top. The 
location of Ho and [N II] lines are not contaminated by strong OH 
lines nor strong telluric absorptions. We achieve a S/N of 37 on 
Ha and of 6.4 on [N II] . 

is con strained to be the theoretical value of lOsterbrockl 
(1989). We fit the line profiles using a weighted \ 2 mini- 
mization procedure which takes into account the greater 
noise level close to atmospheric OH emission. A blank re- 
gion of sky within each reduced datacube was used to de- 
termine the variance spectrum. Since the variance spec- 
trum is dominated by sky emission, the fitted spectra 
were weighted by the inverse of the variance spectrum 
(dominated by sky counts) so that regions of strong sky 
lines do not cause spurious fits. Fortunately, the sky 
emission is weak over the Ha, [Nil] wavelength range for 
A2667 arcA (Figure [2|). We compute the \ 2 statistic for 
the best-fitting Gaussian lines and compare with that of 
a featureless spectrum. We first fit the spectrum for each 
pixel, requiring a minimum of 5a for detection. Where 
this criterion is not met, we average the surrounding 3x3 
spatial pixels to achieve higher signal to noise. No fit is 
made if the 3x3 averaging fails to produce the minimum 
S/N. Figure [2] shows the integrated spectrum from the 
image plane. 

In order to check that our results are insensitive to 
the lensing model uncertainty, we first analyze the data 
on the image plane, and then repeat the analysis on the 
source plane. We find that the results are consistent for 
the image-plane and source-plane analysis. In the follow- 
ing sections, we will show the spectra measured on the 
image-plane with associated physical scales determined 
from the source plane. 

3. GLOBAL PROPERTIES OF THE GIANT ARC 

We determined the total photometry for arc A 
in 7 bands from HST and Spitzer archival data: 
HST ACS/F814, WFPC2/F606W, WFPC2/F450W, 
and Spitzer IRAC 3.6, 4.5, 5.8 and 8.0^m. Chandra data 
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are also available for the cluster, but we find no X-ray 
detection on the giant arc. We used a fixed aperture of 
diameter 4" defined in the ACS/F8f4 band to measure 
the photometry. For the IRAC colors, we measured the 
PSF near arc A using the bright isolated stars in the field 
of view and convolve the ACS/F8f4 aperture with this 
PSF to calculate the aperture flux. Note that the photo- 
metric aperture covers the entire image of arc A whereas 
the OSIRIS FOV does not cover the tip of the southern 
clump C3 (Figure Q] for OSIRIS FOV; C3 is labeled on 
Figure 3). This southern clump only contributes to ~ 
5% to the total photometry, well below the ~10% pho- 
tometric measurement error. The intrinsic photometric 
error is dominated by the lensing magnification uncer- 
tainty (~15%). 

We use the software LE PHARE (jllbert et all 120091 ) to 
determine the stellar mass. LE PHARE is a photometric 
rcdshift and sim ulation package based on po pulation syn- 
thesis models of iBruzual fc Charlotl (120031). We ch oose 
the initial mass fu nction (IMF) by IChabrierl ([2003] ) and 
the lCalzetti et al.l (|2000D attenuation law, with E(B - V) 
ranging from to 2 and an exponentially decreasing 
SFR (SFR oc e"*/ 7 ") with t varying between and 
13 Gyrs. After correcting for the flux magnification 
of fi =14±2.1, the stellar mass from SED fitting is 
Mstar=1.90if;§ x 1O 1O M (lensing uncertainty included), 
assuming the best fit extinction value E(B — V)=0.6, 
and SFRsed = 50± 35 M yr" 1 . T aking E(B - V) gas = 
E(B - V) star / 0.44 (lCalzettil H997h and Rv=3.1. we ob- 
tain A(Ha)=3.46 according to the lCardelli et al.l (|1989f ) 
classical extinction curve. Note that the A(Ha) is ~ 
2.46 higher than the median A(Ha) ( ~ 1.0) at similar 
masses from the local and z ~ 1.7 star fo rming galaxy 
dust extinction vs. stell ar mass relation (jGarn fc Best! 
2010; iSobral et alJ[2012ah . 

The total Ha flux from our OSIRIS spectra is 6.3±0.2 
xlO -16 ergs s _1 cm -2 (befo re correcting for le nsing mag- 
nification). Applying the iKennicuttl JI998) SFR pre- 
scription and after correcting for lensing magnification 
and extinction (A(Ha)=3.46), the Ha star formation rate 
is estimated to be SFR# Q = 49.2±16.5 M Q yr" 1 . Note 
that the error is calculated by assuming a 15% uncer- 
tainty in lensing magnification and 30% error in Ha flux 
measurement. The error estimation excludes the uncer- 
tainty in extinction and aperture flux loss due to the cov- 
erage of the OSIRIS instrument. The physical properties 
of Arc A is summarized in Table 1 . 

4. LENSING RECONSTRUCTION, MORPHOLOGY 
AND KINEMATICS 

4.1. Lens Modeling 

We use the best fi t lensing mod e l par ameters from 
iCovone et all (120061) : iRichard et all dgOlOh and use the 
software Lenstool 5 (jKneib et al.lll993HJullo et al.ll2007l) 
to derive the geometrical transformation between image 
plane coordinates and source plane coordinates. The 
HST image and IFU datacube are reconstructed using 
this mapping and based on the conservation of surface 
brightness. The total magnification is computed by the 
ratio of the sizes (or equivalently, the total fluxes) be- 
tween the image and its source plane reconstruction. As 

5 http : //www . oamp . f r/cosmology/lenstool 



Fig. 3.— Left: HST F814W, F606W, F450W three-color image 
for A2667 arc A. The black contours show the observed Ha flux 
levels of 2, 3, 3.5 in units of 10 — 17 ergs s — 1 cm -2 arcsec - 2 . Right: 
reconstructed source plane morphology from the three-color HST 
image on the left. The physical scale is marked. The distance be- 
tween the center of the red and blue clumps is ~ 600 pc. There are 
three main clumps revealed in both the observed HST and recon- 
structed image: a northern red component (marked as CI), and 
two blue southern clumps (marked as C2, C3). The morphology 
implies either a face-on clumpy disk or a merging system. North 
is up and East is to the right. 

the magnification is not isotropic, the angular size of each 
image is more highly stretched along a specific axis. How- 
ever, within the single image of arc A, the distortions 
among each pixel is identical. Figure [3] shows the recon- 
structed image from the HST F814W, F606W, F450W 
three-color image. The total flux magnification for arc A 
is /i=14.0±2.1, and linear magnifications in two axes are 
jtil=7.0±1.4, and ^2=2.0±0.4. 

4.2. Morphology & Kinematics 

The HST three-color image reveals at least three main 
clumps (labeled as CI, C2, C3) for arc A. These three 
clumps are easily identified in both the image-plane (Fig- 
ure [3j Left) and source- plane (Figure [3l Right) morphol- 
ogy. The northern red clump (CI) could be dominated 
by an old stellar population. Alternatively, CI could 
be heavily reddened by dust. After masking out the 
two blue clumps (C2, C3) on the south, a Sersic pro- 
file fit to the red component on the source plane yields 
a Sersic index of 0.63 ± 0.30, and a half light radius 
of Re = l.l±0.2 kpc. The blue clumps (C2, C3) are 
likely to be active star forming regions. The peak of 
the Ha emission is close to the largest blue clump as 
shown in the Ha contours of Figure[3] One the one hand, 
the asymmetric clumpy structure of A2667arc is similar 
to the asymmetric clumpy feature seen in z ~ 2 disks 
dGenzel et al.ll2006l[2bllHLaw et al ll2009l: Eehnert et al.l 
I2009t iForster Schreiber et alj|2009l |2011[ ). On the other 
hand, the morphology could also be interpreted as a 
merging system. We combine the morphology with the 
velocity structure to further analyze the disk-versus- 
merger origin of the system. 

In Figure HJ we show the two-dimensional (2D) Ha 
intensity, Ha Gaussian fit velocity and velocity disper- 
sion maps. In Figure we extract the one-dimensional 
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FlO. 4. — From left to right: the Ha intensity, Ha Gaussian fit velocities and velocity dispersions presented on the image-plane. The Ha 
flux is in units of 10 -17 ergs s — 1 cm -2 , and the velocity is in unit of km s — 1 . For the middle panel: the black dashed lines indicate the 
"slit" used for extracting ID spectrum in Figure 5; the blue and red boxes are used in the analysis of Section 5.2. These images are shown 
in IFS datacube coordinates. North and East are 37° clockwise from the HST images in Figures 1 and 3. 
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and extract the Ha flux weighted spectrum for each 
bin. There is a small velocity gradient (^75 km s _1 ) 
across a 2 kpc length. Since the shape of the fD 
velocity curve is not clearly of a rotating disk, the 
small velocity gradient could result from either an al- 
most face-on rotating field or a non-rotating system. 
The system however is clearly dominated by disper- 
sion with cr median /Vmax sini~1.3, consistent with the 
large velocity dispersion of the high redshift turbulent 
disks (iLehnert et al.|[2009t iForster Schreiber et al.ll2009t 
I Jones et al.ll2010bft . 

From our morphology and velocity maps, it remains 
unclear whether the A2667 arc is a face-on clumpy disk 
or a merging system, especially considering that the Ha 
lines are not observed across the entire galaxy. Higher 
S /N data with better spatial coverage is required to more 
reliably model the velocity field. 



Fig. 5. — One-dimensional velocity extracted from the slit po- 
sition shown Figure 4. The five data bins from left to right cor- 
respond to the direction from box 5 to box 1 in Figure 4. The 
physical scale is calculated on the source-plane. Note that the er- 
rors bars on each data point are the bootstrapped median errors 
for all the pixels in the bins, weighted by the Ha intensity. The 
high velocity dispersion is consistent with those of turbulent high 
redshift galaxy disks. 

(ID) velocities and dispersions from a slit oriented along 
the direction of the blue and red components (black dot- 
ted lines in Figure 4, middle panel). We divide the 
spaxels into five bins equally along the slit direction 



5. SPATIALLY RESOLVED LINE RATIOS 

5.1. [Nil] \6583/Ha as an indicator for metallicity 
and nonphotoionized emission 

The [Nil] A6583/Ha ratio is often used as an indicator 
for nebular metallicity at high redshift (up to z ~ 2.5). 
The adjacency of the two lines in wavelength means the 
line ratio is insensitive to redde ning correction a n d flux 
calibration. The calibration of iPettini fc Pagell ()2004f ) 
(PP04) is commonly used to convert [N II] /Ha to metal- 
licity. Using the oxygen abundance of local H n regions, 
PP04 measured metallicities from the Te (electron tern- 
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Fig. 6. — The [Nil] /Ha ratio variation in boxed regions along the major axis. The top 6 panels show the observed image-plane spectra 
in vicinity of the [N II] and Ha for the 6 boxed regions in Figure 4. The spectra are the Ha flux weighted median of the pixels in each box. 
The x axes are wavelengths in A, and the y axes are flux density are in 0.001xl0 — 17 ergs s -1 cm -1 A -1 . The [Nil] and Ha line positions 
are marked as vertical dashed lines. The bottom panel shows the [N II] /Ha ratio as a function of the intrinsic physical scale (after lensing 
correction) of the galaxy. The red horizontal line indicates the line ratio corresponding to the solar metallicity value. Ratios above solar 
can not be directly related to metallicity. We also notice a significant enhancement of [Nil] /Ha in Box B, a blueshifted region with respect 
to the galactic center. 



perature) method and determined the correlation be- 
tween metallicity and the N2 (= log([N n] A6583/Ha)) 
values as follows: 



12 + log(0/H) = 8.90 + 0.57 x N2 



(!) 



Note that at high abundance, PP04 used the strong 
line metallicities rather than the Te metallicities because 
of the weakness of the Te sensitive [O in] A4363 line. 
Strictly speaking the PP04 relation only holds true for 
the range of -1.8 < N2 < -0.8 (or 0.0158 < [Nil] /Ha 
< 0.158). At extremely low N2 values, the calibration 



of N2 is difficult due to the weakness of the [N ii] line. 
At N2 > —0.8, the relation is nonlinear. The reason for 
this is, first, as the metallicity becomes larger and Te 
decreases, the production of secondary nitrogen begins 
to dominate and maintain a high value of N2. At even 
higher metallicity, the rise in metallicity cools the nebula, 
lowering Te and allowing less collisio nal excitations of the 
[Nil] lines (|Kewlev fc Dopital l2002j ) . We note that the 
[N II] /Ha ratio is particularly sensitive to non-thermal 
excitation. The presence of shock excitation or an AGN 
will increase the [N ii] /Ha ratio and cause abundances 
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determined to be artificially high ( Baldwin et al.l 119811 : 
iKewlev fc Dopital l2002t iRich et alJ [2010[). For metallic- 
ity studies, it is crucial to check if the N2 values are 
within the regime of pure star formation. 

On the other hand, if N2 surpasses the linear metallic- 
ity diagnostic range, the high N2 ratios may be a smoking 
gun for the presence of non-ph otoionization mechanisms 
such as shocks. For example, IVeilleux fc Rupkel (2002) 
studied the [N n] /Ha distribution in the local super- wind 
galaxy M82. They found that for [Nil] /Ha > 0.5, addi- 
tional excitation mechanisms such as shocks are needed 
to explain the enhanced [Nil] /Ha ratio. Using the most 
recent sho ck and photqioniz ation models with IFU ob- 
servations, IRich et all ([20111 ) showed that for the local 
galaxy IC 1623, shocks contribute more than 40% of the 
line ratio of [Nil] /Ha. 

5.2. [Nu] /Ha in boxed regions 

To determine the spatial variation (if any) of metallic- 
ity or ionization sources, we analyze the spatial distribu- 
tion of [Nil] /Ha ratios. In order to enhance the S/N of 
the N2 line ratio map, we bin the pixels in five rectan- 
gular regions along the major axis and extract spectra 
for each boxed region (Figure 0}. The widths of the re- 
gions are chosen in such a way that each region contains 
sufficient pixels to extract a reliable spectrum (i.e., S/N 
> 3 for the [Nil] line). We derive the spectrum for each 
region using the Ha flux weighted median of all pixels in 
that region. Wc fit the [N n] and Ha emission lines and 
obtain the errors on [N n] / Ha by propagating the flux 
errors of [N n] and Ha. 

The peak of the Ha emission corresponds to the cen- 
ter of Box 3 in Figure |U We do not find significantly 
enhanced [N n] /Ha ratios in the central pixels and rule 
out AGN contamination of this target. 

Our [N n] /Ha emission-line ratios are shown in Fig- 
ure El The [Nil] /Ha ratios of Box 3 and 4 are 
0.303±0.034 and 0.312±0.033 respectively, close to so- 
lar values. The [N n] /Ha ratios for Box 1 and 2 arc 
0.472±0.046 and 0.384±0.03 respectively, significantly 
higher than those of Box 3 and Box 4. We can only 
derive an upper limit of [Nil] /Ha (<0.233) for Box 5. 
The high [N n] /Ha values in the non-central regions of 
the galaxy indicate either extremely high metallicity or 
the presence of non-photoionization excitation. 

We further analyze individual pixels in Box 1 and 2 
and find that the [N n] /Ha ratio is the highest in a sub- 
region highlighted as Box B in Figure 2J Very inter- 
estingly, Box B which has the highest [N n] /Ha ratio 
([Nil] /Ha =0.499±0.05) has also the "bluest" velocity 
(Av ~ 150 kms" 1 ) with respect to the center of Box 
3. Unfortunately, the data quality is not sufficient to 
reveal any significant variations in the velocity disper- 
sion map. As discussed in the previous section, large 
[Nil] /Ha ratios may indicate non-photoionized emission 
and therefore we do not convert this ratio to metallic- 
ity. We note that the consequences of deriving metal- 
lically using [N n] /Ha ratios in this case will result in 
the conclusion of a flat or inverted metallicity gradient 
(Appendix). In the next section, we discuss the possible 
origins of these elevated [Nil] /Ha ratios. 

6. DISCUSSION 
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Fig. 7. — The starburst-shock mixing grids on the [Om]/H/3 
versus [N II] /Ha diagram. The dashed curve labele d as "KeOl" is 
the the oretical "maximum starburst line" derived bv lKewIev et all 
H2001al) as an upper envelope for st ar-forming galaxie s. The dashed 
curve labeled as "Ka03" is the Kauffmann ct al. (2003) semi- 
empirical boundary for star-forming galaxies. The pure HII region 
models are shown as the blue solid line and the pure shock models 
are shown as the red solid line. The colored dots in between form 
a mixing sequence from pure HII region to pure shock excitation. 
Different colors represent different shock velocities and the relative 
contribution of shock is shown in different symbols. Solar abun- 
dance is adopted. The [N II] /Ha ratio for region B in Figure [6] is 
highlighted as the vertical red dashed line. 



The standard BPT diagrams prove to be a pow- 
erful tool for investigat i ng th e ionization sources in 
galax ies (Baldwin et al.l [19811: Vcillcux & Oste rbrockl 

KT 



19871: IKewlev et all I2001at iKauffmann eTaTl [2003; 
Kewlev et al.ll2006l) . Employing IFS data on the BPT di- 



agram has revolutionized our und erstanding of the power 
sources across entire galaxies (e.g . .lRich et ail2 010. 2011, 
|20I3 ISharp fc Bland-Hawthornl |2010j). The [Nil] /Ha 
versus [O in] /H/3 diagram is useful for distinguishing 
pure star-forming regions from non-photoionized re- 
gions. The [N 11] /Ha ratio increases linearly with 
nebular metallicity until high metallici ties where the 
log([Nn1 /Ha) ratio saturates at -0.5 (iDenicolq et all 
120021 IKewlev fc Dopital 120021 : iPettini fc Pagell 12004). 
Any non-thermal contribution will shift the [N 11] /Ha 
ratio towards higher values than this s aturat ion level 
(|Kewlev et al.1 [2006L IRich et al.l [2fTinL l20ll . The 
log([Nn] /Ha) ratio of region B is 0.2 dex larger than 
the saturation level of —0.5, suggestive of shock ioniza- 
tion. 

We apply a mixed grid of photoionization and shock 
models to quantify the potential relative contribution 
of shock and H 11 region ionization to the [N 11] /Ha 
ratio. The Hn region models are taken from the 
most recent Starburst99(SB99) /Mappings III model 
suite ([Leitherer et al.lll999fc [Sutherland et al.lll993l) . The 
m odel setup and i nput p arameters are similar to those 
of iLevesque et all (|2010| ) except that we use a slightly 
higher mass cut of 120 M Q for the Salpeter initial mass 
function (Salpeter 1955). We use a continuous star for- 



mation history with age 5.0 Myr and n e =100 cm~ 3 . So- 
lar metallicity was used with a range of ionization param- 
eters log U (—3.5 to —1.9). We then mix the Hn region 
model with the slow- vel o city shock model as described in 
ISutherland. et all (|1993l ): lFarage et all (|2010f ): lRich et al.l 
(|2010L [201 ll ) . We assume a 100% pre-ionization, with a 
slow shock velocity ranging from 100 to 200 km s _1 , in 
increments 20 km s - 1 . These shock velocities are typi- 
cal of galac tic winds (Shopb ell fc Bland-Hawthornlll998l : 
iRich et al.ll2O10L 120111: iSharp ~B land-H awthornll2010h ' 

Our shock and star burst models are shown in Figure [3 
We show mixing lines to represent the relative contribu- 
tions of starbursts and shocks from to 100%. Shocks 
are able to produce higher [Nil] /Ha ratios than pure 
starburst models for a given metallicity. According to 
our models, a 20% shock contribution could increase the 
log([Nn]/Ha) value by 0.2 dex. The region B of the 
A2667arc has log([N n]/Ha)=— 0.3, lying well within the 
region covered by the shock models. 

Shocks induced by galactic winds/outflows as- 
sociated with either stellar winds, supernovae or 
nuclear activi ty are commonly observed in local 
galaxies (e . g.. iVeilleux et al.l 12005); iRupke et al.l 120051: 
IRich et all 120101: iSoto fc Martini 120m The ubiq- 
uity of winds is even more prominent at high 
redshift gi ven their violent s tar formation activi- 
ties (e.g. iPettini et al.l 120011: iSwinbank et ail 120051: 
Sato et al.1 120091: IWeiner et al.l l2009t iMartin fc Bouchel 



20091: IMartin et al.l l2012HRubin et al.1 l2010bL I20TTF 
Alexander et~aT] 120101: fSteidel et al.l 12010ft . In the cold 
flow sc enario, the turbule nt disks are associated with 
shocks ()Genzel et al.|[201ll ). The coincidence of region 
B being blue-shifted with respect to the center strength- 
ens the argument that it is an outflow shocked region. 
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Fig. 8. — The co-added spectrum of Boxes 1, 2, 4, 5 (i.e., non- 
nuclear regions) in Figure 6. The x axis is wavelengths in A, and 
the y axes are flux density are in units of 10 -17 ergs s — 1 cm -1 A - 1 . 
The top panel fits a single component to the Ha and [N ll] lines 
(blue solid curve). The bottom panel fits a narrow component 
(green curve) to the Ha, [Nil] lines and adds a broad component 
(red curve) to the Ha line; the sum of the two-component fit is 
shown (blue curves). The broad Ha component on the bottom 
panel (red line) is ~ 60 km s — 1 blue-shifted from the narrow-line 
component. 



Indeed, the outflow scenario is further supported by 
the detection of a broad component in Ha. We detected 
a broad component in the co- added spectrum of boxes 1, 
2, 4, 5. In Figure [8] we show that by including a broad 
component in Ha, the fitting to the co-add spectrum is 
improved by 3<r in a \ 2 sense. The broad component 
indicates a broadening wing of Av(FWHM) ~ 380 ± 130 
km s _1 , with a slight blueshift of ~ 60 km s _1 with 
respect to the narrow line component. The contribution 
of the broad wing to the total Ha flux of the co-added 
spectrum is f broad ~ 40%. 

Because of the unknown orientation of the A2667arc 
system, the outflow velocity is difficult to estimate. We 
thus compare the outflow properties of the A2667arc with 
other studies using the observed linewidth Av(FWHM). 
For example, NGC 1569, a local dwarf starburst galaxy, 
has a wind-e xcited broad component of Av(FWHM) < 
300 km s" 1 (jWestmoauette et al.ll2007t ). Arp 220, a fa- 
mous local IR-luminous starburst galaxy, has a nuclear 
wind broadening of Av(FWHM) ~ 300 - 800 km s _1 
(|Arribas et al.l |2001|). For a l arge sample of 39 local 
ULRGs studied in ISoto et al.l ((2012), the median line- 
width for the shocked components is Av(FWHM) ~ 
339 km s _1 (after converting their a to FWHM using 
FWHM= 2 3548 xGaussiana). The median SFR of the 
Soto et al.l (f20ll) ULIRGs sample is ~ 73 M Q yr" 1 . At 
very high redshift (z > 2), the UV-selected massive SINS 
;alaxies are found to have At> (FWHM) ~ 500 km s _1 
Shapiro et all 120091: iGenzel et ail l2011h . with a me- 
dian SFR of ~ 72 M yr -1 . Roughly, the broaden- 
ing linewidth (Av(FWHM) ~ 380 ± 130 km s" 1 ) and 
SFR (49.2±16.5 M Q yr" 1 ) of A2667arc is most similar 
to the local ULIRGs sample of ISoto et all (|2012h . How- 
ever, the similarity between the linewidths and SFR of 
the A2667arc and the z ~ 2 SINS galaxies cannot be 
ruled out within 3cr of the Av (FWHM) measurement 
uncertainty. 

If, on the other hand, the high extra- nuclear [N n] /Ha 
ratio regions in the A2667arc are caused by extremely 
high metallicity, we lack theories to explain how outer re- 
gions with super-solar metallicity would evolve spatially 
as a function of time. It is unclear how the A2667arc has 
gone through such a quick enrichment. Throughout the 
study we use the peak of Ha as the "galactic center", 
which does not necessarily concur with the true galactic 
center. However, the choice of center would not change 
our main conclusion that there is a shocked region blue- 
shifted from the peak of Ha emission. 

Finally, the possibility of shocks will cause an overesti- 
mation of the "true" SFR. The Ha flux in Regio ns 1 and 
2 is -20% of the total Ha flux. IRich et al.l (|2010l ) showed 
that for a slow shock of ~ 100 km s" 1 , 50% to 70% of the 
Ha luminosity can be ascribed to shocks. Hence, ~10% 
of the total Ha flux of A2667arc may result from shocks. 

7. SUMMARY 

We present the OSIRIS IFU data analysis of the spa- 
tially resolved giant arc at z = 1.03 in Abell 2667. With 
a flux magnification of 14.0±2 in flux and — 2 x 7 in 
area, we achieved a spatial resolution of 115—405 pc on 
the source plane. We find the A2667arc is clumpy in 
rest-frame UV and optical morphology, and has a high 
velocity dispersion V ro t smi/a — 0.73, similar to the 
turbulent star forming disks at z— 2. The current mor- 
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phology and velocity map is not sufficient to distinguish 
whether the A2667 arc is a face-on clumpy disk or a 
merger. We find that [Nil] and Ha ratios first rise and 
then fall when moving out from the galactic center de- 
fined as the peak of Ha flux. The central 350 pc of the 
Ha nucleus region is dominated by star-formation with 
metallicity 12 + log(0/H) =8.57±0.03 (~ 0.76 Z Q ). The 
[N n] /Ha ratios in the 350 — 900 pc annulus show a 
significantly enhanced value (~ 0.1 dex higher than the 
nucleus). The highest [Nil] /Ha ratio region (~ 0.2 dex 
higher) corresponds to a ^100 kins" 1 blue-shifted ve- 
locity clump with respect to the nucleus. Applying the 
most recent slow-shock models, we propose that this ele- 
vated [N n]/Ha ratio is caused by shock excitation from a 
galactic wind. Our results indicate that extreme caution 
must be used when interpreting line ratios as metallicity 
gradients at high redshift. Shocked extra- nuclear regions 



can mimic flat or inverted metallicity gradients. 
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Table 1. 

Derived Global Physical Properties for Arc A 



Re 
kpc 


E(B - V) 
SED fit 


SFR Hq! 

Meyr- 1 


Mstai(SED) 
10 10 M 


a 

krns" 


Vmai sin i 
krns' 1 


12 + log(0/H) 
Integrated 


flux, axisl, axis2 


1.1±0.2 


0.6 


49.2±16.5 


1.90^-0 


65±6 


48±12 


8.606±0.051 


14±2.1, 2.0±0.4,7.0±1.4 



Note. — Derived intrinsic global properties for A2667 arc A. Col. 1: The half light radius from fitting a Sersic profile to the red HST 
component. Col. 2: Extinction derived from SED fitting. Col. 3: SFR from Ha flux, corrected for lensing magnification and extinction. 
The error comes from the uncertainty in magnification and flux calibration. Col. 4: Stellar mass from SED fitting. Col. 5: Median velocity 
dispersion, weighted by Ha flux. Col. 6: V m ax sini see the text for methodology. Col. 7: Metallicity from the integrated spectrum. Col. 
8: Flux magnification and size magnification in the major and minor axes. 
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APPENDIX 
[Nil] /Ha AS A FUNCTION OF RADIUS 

In metallicity gradient studies, ave raging spectra in annular regions is a common exercise (|Moustakas et all [2010; 
iQuevrel et al.ll201ll : lYuan et aUl2011l) . To relate to these studies, we extract spectra in 3 annuli centered on the Ha 
peak (Figure IHJ). Each spectrum is calculated as the Ha flux weighted median of all pixels within the annulus. From 
Section 15.21 we see that the line ratio map is not symmetric with respect to the center. Averaging spectra in annuli 
smooths out the lateral spatial variance, resulting in noisier [Nil] lines. The choice of the annuli is such that the 
integrated spectra gives S/N > 5 for the [Nil] lines. The physical lengths of the three outer radii of the annuli in the 
source plane are 0.35, 0.93, and 1.31 kpc respectively. The [Nil] line is robustly detected at S/N > 5 for the inner two 
annuli and is a 3<7 detection for the outer annulus. 

The [Nil] /Ha ratios are larger in the circum-nuclear region (350—930 pc) than in the nucleus (<350pc). Small 
[Nil] /Ha ratios are seen in the outer annulus (930—1.3 kpc). If the [Nil] /Ha ratios arise purely from Hn regions, 
these [Nil] /Ha ratios could t ransla t e into a flat or i nvert ed metallicity gradient similar to those in several high redshift 
studies such as ICresci et~aLl ()2010D ; IQuevrel et al.l (|2011l) . 

However, in A2667arc, these enhanc ed circum-nuclear [N ii] / Ha ratios are remarkably similar to the LINER-like 
luminous infrared galaxies reported in I Rich et al.l (|2010l . 1201 ll ). Rich et al. observed that enhanced extra- nuclear 
[Nil] /Ha ratios in late-stage galaxy mergers are driven by slow shock emission (100—200 km s _1 ) associated with 
galactic-wide winds or gas collisions due to mergers. Rich et al. (2010, 2011) showed that both velocity dispersions 
and the strong line ratios independently indicate shock excitation in the extra-nuclear regions. Slow shock models 
confirm this conclusion, yielding remarkable agreement with the enhanced emission-like ratios. 
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Fig. 9. — The [Nil] /Ho ratio as a function of the "galactoccntric" distance. The galactic center is defined as the peak of Ha. The 
top left panel shows the 3 annuli relative to the H« intensity map. The top right panel shows the 3 annuli relative to the velocity map. 
Note that these images are given in IFS datacube coordinates, thus North and East is 37° clockwise from the up and left position. The 
following 3 panels show the observed image-plane spectra averaged in the 3 annuli. The wavelength is in units of A, and flux density is in 
units of 10 — 17 ergs s — 1 cm -1 A" 1 . The line positions are marked on the plot . The last panel shows the [Nil] /Ha ratio as a function of 
the intrinsic "galactocentric" distance (after lensing correction). The red horizontal line indicates the line ratio corresponding to the solar 
[N II] /Ha ratio. Ratios above solar may be contaminated by non-photoionizcd emission and can not be directly related to metallicity. 



